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Abstract 

The propagator in the instanton background in the (— X(f) ) scalar model in 
four dimensions is studied. Leading and sub-leading terms of its asymptotics 
for large momenta and its on-shell double residue are calculated analytically. 
These results are applied to the analysis of the initial-state and initial-final-state 
corrections and the calculation of the next-to-leading (propagator) correction to 
the exponent of the cross section of instanton induced multiparticle scattering 
processes. 

1 Introduction 

There is a number of interesting physical effects induced by instanton solutions. They 
appear in theories which possess a non-trivial structure of vacua. The most promi- 
nent example is the electroweak theory with an infinite number of vacuum states la- 
belled by the Chern-Simons number The instanton solutions describe transitions 
with baryon number violation between the vacua. Another example of the instan- 
ton induced process is the decay of a metastable (false) vacuum due to underbarrier 
tunelling from a false vacuum to the true one f2j. The third example is a shadow 
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process This is a non-perturbative process in which both the initial and the final 
state are in the false vacuum. Apart from standard perturbative contributions, the 
processes which start and end in the false vacuum acquire additional contributions 
due to the underbarrier tunelling of the system to another vacuum and its return 
to the initial one. This transition is obviously induced by an instanton solution and 
goes through the intermediate state containing a bubble of the true vacuum. In this 
article we study the cross section of shadow processes. 

Much work has been done to study of the instanton induced transitions, and quite 
effective techniques for the calculation of the probabilities of such transitions have 
been developed (see Refs. 0, [5| for a review). Consider for simplicity a scalar theory 
with the field and the action S ((/)). Let us discuss a process (2 — > any) with two 
initial particles of the total energy E. The total cross section of this process is given 
by the sum over the partial cross sections, 

n 

and each term is calculated from the amplitude A 2 ^ n (pi,P2] Qi, ■ ■ ■ , Qn) in a standard 
way. The amplitude is obtained by applying the Lemann-Symanzyk-Zimmermann 
reduction formulas to the (n + 2)-point Green function, 

G n+2 ( Xl , . . . , x n+2 ) = J D0(x)e- sw 0(^i) • • • <K*n +2 ). (1) 

If the model possesses an instanton solution 4>i ns t(x), then besides the standard pertur- 
bative contribution to the Green function, corresponding to the expansion of Eq. ([!]) 
around the trivial solution = 0, there exists a contribution due to the instanton 
sector. It is the contribution we will focus on in this paper. It can be calculated by 
performing the expansion (j){x) = <pi ns t{x) + t](x) and integrating over the fluctuations 
r](x) around the instanton: 

G n+2 (x 1 , x n+2 ) = J Dr](x)e~ Smst exp j-~ J dxr](x)D x r](x) - ... 

X [(pinst(Xi) + T}(Xi)] . . . [(t>inst(x n+2 ) + 7?0„ +2 )] • (2) 

Here Si ns t is the instanton action and D x is the second order differential operator 
determined by the quadratic in r](x) part of the action. The dots in the exponent 
stand for cubic and higher order terms in rj(x). In general the instanton solution 
is parametrized by a number of continuos parameters denoted here by (a- These 
parameters correspond to symmetries of the model. Thus, in fact, there is an infinite 
family of instanton configurations which have to be taken into account. Due to this 
fact the operator D x possesses zero modes associated with these symmetries, and is 
not invertible. Functional integration over the fluctuations r](x) is perfomed according 
to a known procedure, and, as a result, instanton contribution (Q) to the Green 
function G n+2 (xi, ■ ■ ■ iX n+2 ) includes the integral over the parameters (a (see ref. [|J 
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for details). In the leading semiclassical approximation the result can be schematically 
written as 



G n+2 (x 1 , x n+2 ) ~ / d( A e Smat (j)inst(xi; C) • • • <t>inst{xn+2] C) + corrections. (3) 



Using Eq. (|3|) one can evaluate the partial cross-sections cr 2 ^ n (E), and after per- 
forming the summation over n one obtains the leading-order semiclassical expression 
for the total cross-section: 

a 2 (E) ~ exp {1 (F^(E/E sph ) + ...)+ 0(1)} , (4) 

where A is the coupling constant in the model and E sp h is the energy of the sphaleron 
configuration which characterizes the height of the barrier separating the vacua. The 
superscript (L) stands for " leading" . On dimensional grounds it can be shown that 

Tfl 

E sph = k—, (5) 

where k is a numerical constant depending on the model. It is easy to see that for a 
shadow process 

F^ L \0) = 2XS inst . 

The dots in the exponent of the cross section in Eq. (f|) stand for A terms. In the 
weak coupling regime, A — > 0, contributions described by this function dominate and 
the 0(1) term in Eq. (|4]) gives an exponentially small correction. 

The first calculations of instanton induced processes were done for the electroweak 
theory in Refs. 0. There the role of A is played by g 2 , the square of the gauge 
coupling constant and Si nst = 8n 2 /g 2 . The leading order contribution to the function 
F(E/E sp h) is equal to 
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where the parameter Eq = \^6M W / aw is of the order of the sphaleron energy E sph w 
lOTeV0,| 

In the scalar (— A0 4 )-theory, which we are going to consider in this paper, the 
leading contribution is equal to 

where e = E/E sph and E sph is given by Eq. (|5]) with k = 113.4 |fj, 

An effective method of calculation of the corrections to the function F(E/E sp h) 
was proposed and developed in Refs. M. 
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From the above examples it is clear that at energies E <C E sp h the cross section of 
the instanton induced process is exponentially supressed. On the other hand there are 
reasons to expect that as the total energy of the initial states grows, the probability of 
the underbarrier transition, induced by the instanton, grows as well and the supression 
becomes weaker. For high energies the function F(E/E sp h) may be close enough to 
zero and the instanton induced processes may become observable. This opens a 
possibility of new interesting physics within the Standard Model at future colliders 



T0|| • However, in Ref . 1 1 LJ arguments were presented which suggest that perhaps the 
function F(E/E sp h) never vanishes and, hence, at all energies some supression always 
remains. 

The considerations above are based on the assumption that the complete expres- 
sion for the total cross section, not just its leading order contribution, can be presented 
in the exponential form, i.e. 

a 2 (E) ~ e^W+oU (7) 

The leading order term in F has been just discussed. The next-to-leading term is a 
propagator correction, for after integration over i](x) in Eq. (||D it includes contribu- 
tions from the propagator in the instanton background. We will call this propagator 
the instanton propagator for shortness. It is, of course, equal to D^ 1 restricted to 
the space of fields orthogonal to the zero modes. In fact, as we will see, calculation 
of the next-to-leading correction requires knowledge of the exact expression for the 
residue of the instanton propagator. It turns out that in the (— A0 4 )-theory an exact 
expression for the instanton propagator can be obtained. Derivation of the exact for- 
mula for the residue of the instanton propagator, as well as calculation and discussion 
of the propagator correction to the function F(e) in the scalar model is one of the 
purposes of the present paper. 

An important issue in the theory of instanton induced processes is the validity 
of formula (0). The difficulty in proving it is due to the initial state corrections and 
initial-final state corrections. In general, it is not easy to show that they exponentiate, 
so that the whole result can be represented in the form given by Eq. ([F]). In the 
electroweak theory a proof, based on the properties of the propagator in the instanton 



background was given in Ref. |]12|| . We apply the arguments of Ref. |12| in the 
(— A0 4 )-theory making use of the explicit expression for the propagator. 

On the other hand, it was shown (see Refs. |13|, [L4|, |I5|1 ) that for the (N — > any) 
process, where N ~ 1/A and, therefore, is parametrically large for small A, the total 
transition probability in the semiclassical approximation is given by 

<r N (E) ~ e&^+oW. (8) 

Here we introduced the parameter v = N/N sp h, where 

N sph = j (9) 
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is the characteristic number of particles contained in the sphaleron. In Refs. |13|, [14] 
it was conjectured that in the leading semiclassical approximation the two-particle 
cross section can be calculated from the following formula: 

lim Alna 2 = \imF(e,u) , (10) 

A^O v— >0 

provided the limit v — > exists. The problem with this conjecture is that the function 
F(e, v) is known to contain contributions singular in v. In particular, such contri- 
butions appear in the course of the calculation of the propagator correction. The 
conjecture of Refs. fl3| , [T4}| basically claims that all such terms cancel each other 
in the final answer. Its validity, of course, means that the semiclassical form of the 
two-particle cross section is indeed given by Eq. (^) with F(e) = F(e, 0). Verification 



of conjecture (|To|) in the next-to-leading order is another purpose of this paper. Note 



that different arguments in favor of this conjecture were given in Refs. [IIB, [T? 



The plan of the article is the following. We start Sect. |2] with a review of basic 
facts about the perturbation theory around the instanton following mainly the results 
of Refs. [0], Hl4| . We also discuss the general structure of the leading and next- 



to- leading corrections to the function F(e, v) and saddle point equations defining 
these corrections in the semiclassical approximation. Finally, using a formula for the 
leading asymptotics of the instanton propagator, obtained in Ref. [113, we explain the 



appearance of terms singular in v as v — > and prove that they cancel each other. In 
Sect. H] we describe the model and review the derivation of the exact explicit expression 



for the instanton propagator obtained in this model some time ago in Ref. |T9 |. In 
Sect. H we calculate three leading terms of the large- s asymptotics of the Fourier 
transform of the propagator and discuss the implementation of Mueller's idea in the 
scalar model. In Sect. |5| we obtain the exact expression for the double on-mass-shell 
residue of the instanton propagator. In Sect. |6| we discuss the leading and propagator 
corrections to the exponential of the cross section in the limit v — > 0. An analytical 
expression for them is obtained. We also discuss the terms in v for this concrete 
example. In Sect. |^ the propagator correction to the function F(e, v) for a wide range 
of e and v is calculated numerically and analyzed in detail. In Sect. § we discuss the 
consistency of the procedure used for the calculation of the propagator correction. 
Sect. H contains some discussion of the results and concluding remarks. 



2 Perturbation theory around the instanton 

Within the formalism of coherent states the cross section ([§D can be calculated 
using the formula 

a N (E) = T^pY. I< b\SV E V N \a >| 2 , (11) 

J a,b 

where S is the S'-matrix in the one-instanton sector and Ve and Vn are the projectors 
onto the subspaces of fixed energy E and multiplicity N, respectively. The coherent 
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states are characterized by complex variables a k , &k, where k is the spatial momentum 
and the summation over coherent states in Eq. (|ll]) stands for the functional integral 

y\ . . = J Da^Da^Db^Db^ exp (— J dka k a k — / dkb^b^ 

a,b 



Using results of Refs. m, 14 1 for the matrix elements of the S'-matrix and the 



projectors Ve, Vn in the coherent state representation we find that 

a N (E) = J Da^DalDb^Dbldpdp d A x d A ^d6 exp [-a* a - b*b - iP£ - iN9] 

r p (b* k e lkx °, aice -^o+ik^ + r * (6k; a * } ] ^ (12) 



x exp 



where kx = «Jk(^o)o — kx , w k = y/k 2 + m 2 . Among the parameters of the family of 
instanton solutions we explicitly indicate the position of the center of the instanton 
denoted here by xq, whereas p stands for the rest of the parameters. The parameter 
Xq corresponds to the translational symmetry of the theory. The term (— 1P£), = 
(E, P) is the total 4-momentum of the initial state, comes from the matrix element of 
the projector Ve- Correspondingly, the term (—iNO) in the exponent appears due to 
the projector Vn- See Refs. 0, for the details. The "effective action" r p (& k ,a k ) 
to the leading and next-to-leading orders has the following form: 

r P (& k ,ak) = -S ins t{p) + ri(6 k ,a k ) + r 2 (& k ,a k ) + ■ ■ ■ , (13) 

ri(6 k ,a k ) = aR a + b*R b + ab*, (14) 
T 2 (& k , «k) = ~aR aa a + aR ab b* + -b*R bb b*. (15) 

The dots in ( O ) stand for higher order terms in a k , 6 k . The integration over the 
spatial momentum is implicit, so that aR a = / dka^R a {k), etc. 

The functions R a (k) and R b (k) in Eq. ( |14|) are proportional to the residue of the 
Fourier transform of the instanton field continued to the Minkowski values of k and 
taken on the mass shell. Their definitions are the following: 

Ra(k) = j^^-^=(k 2 + m 2 )(p inst (ko,k;0,p)\ ko=iWk , (16) 

RbO*) = (2^)3/2 v^^ ^ 2 + m2 )^n S t(k , -k;0, p)\ ko= _ iWk . (17) 
The Fourier transform of the instanton field is defined in the standard way: 

4>inst(k; x 0l p) = / d 4 x4>i ns t(x; x ,p)e tkx , (18) 



where we explicitly indicated the dependence of the instanton solution on xq and p. 
The function R a (k) corresponds to the initial state, whereas R b (k) corresponds to 
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the final state. Let us denote the residue of the Fourier transform of the instanton at 
x as I(p): 

I(p) = (k 2 + m 2 )(j) inst (k;0, p)\ k 2 = _ m 2. (19) 
It is clear that for a scalar theory I(p) is independent of k and 



The functions _R aa , i? a t and _R;,& in ( p~5| ) are defined through the double on-mass- 
shell residue of the Fourier transform of the propagator in the instanton background: 



i? aa (k,q) = — — - = {k 2 + m 2 )(q 2 + m 2 )G(k , k; g , q; 0, p) | kQ= 



(2n 




1 


1 


(2tt 




1 


1 


(2tt; 





<j =mi q 



^ab(k,q) = /0 _,, - - ;(fc 2 + m 2 )(g 2 + m 2 )G(fc ,k;go,-q;0,p)| kn ,(21) 

go=— iWq 



R bb (k,q) = —(k +m)(q + m )G(k Q , -k; q Q , -q; 0, p) \ ko 



qo=— ivlq 



The instanton propagator G(x, y; Xo, p) is determined by the operator D x of quadratic 
fluctuations around the instanton (see Eq. (|2|)). Calculation of the propagator and 
of its residue in the scalar model will be the main subject of discussion in Sects. 0-|5[ 
The Fourier transform of G(x, y; Xq, p) is defined in the standard way: 

G(k, q; x , p) = J d A xd 4 ye lkx+ ^G(x, y; x , p). (22) 

Again, R aa corresponds to the propagator connecting particles of the initial state, 
R ab connects a particle of the initial state with a particle of the final states and R bb 
connects two particles of the final state. It is easy to see that the dependence of the 
double on-mass-shell residue of the propagator on the momenta k, q is of the form 
(k, q). In a scalar model the momentum dependence of the Fourier transform of the 
instanton propagator is simple. For example, one can choose k 2 , q 2 and s = (k + q) 2 
as variables. Then the double on-mass-shell residue of the propagator for all three 
formulas fl2~l~|) can be expressed in terms of a single function P(s; p), and one gets 

R *^ ) = J^2^ P M k '^- < 23) 

where # = aa, ab, bb, the s-variables are defined by 

s aa (k, q) = s bb (k, q) = -2m 2 - 2(^ k ^ q - kq), (24) 
s ab (k, q) = -2m 2 + 2(^ k u; q - kq). (25) 

in accordance with definitions (pl|) , and the function P is given, for example, by 

P(s aa (k,q);p) = {k 2 + m 2 ){q 2 + m 2 )G{k,q;0,p)\ ko ^. (26) 

qo=iWq 
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The integral over the variables a^, a£, 6k and 6 k in the leading and next-to- leading 
approxmations is Gaussian and can be readily performed. Indeed, we can write the 
part of fll^D , including these variables, as 



Dipy_ exp 



(27) 



where ip£ = (a^, a£, b^, and the matrix M and the source J can be read from 
pi), PD - (pi). In particular, 



M = M (1) + M (2) + . . . 



(2* 



where M^, corresponding to the leading contribution, includes the terms that do 
not contain the propagator residue, whereas contains propagator contributions 
of Eq. flrSp. As we will see later, in the scalar model ~ (mp) 2 , where p is the 
size of the instanton. For the range of energies that will be studied in Sect. |6| and 
Sect. [7| {mp) is a small parameter. Then ~ (mp) 4 is indeed a small correction 
comparing to the leading term. In the leading approximation integration over ip gives 



Z ~ exp 

where ip is the saddle point value 



exp 



-J t (m^ 



)r 1 J 



j. 



(29) 



calculated in the leading order approximation. Then the cross section in the leading 
approximation is equal to p4 



where 



dpdp'd^xod^^dO exp 



-Sinstip) ~ Sinst(p') + -rW^(x ,p,p',C, 



±W® (x , p, p', e, 0) = -iPi - iNO + l -J T (M«) 1 J 



(30) 



(31) 



T 



XT 



-i 



Here we introduced the variable 



- i P C -iN9 + Rt T -^R b + 1 R: T -^ 



Ra 



7 X 



7 = exp(tu) 

and the functions 

X(k, q) = 6{k - q)e™ k5 , T(k, q) = 5(k - q)e™ ka;o , 



(32) 
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which are treated as matrices in the momentum space, so that the integrations in Eq. 
(]31|) are viewed as matrix multiplications. 

In the next-to-leading approximation the cross section is given by 



<Jn{E) 



dpdpd xod £d9ex.p 



-s mst ( P ) - s inst ( P ') + jWW(x Q ,p,f/,z,e) 



(33) 



where the contribution is obtained by evaluation of the term 



(see (p7D , fl28|) ) at tj) = ijj given by the leading order value (p9|) . This is sufficient for 
the accuracy with which the next-to-leading order is calculated. One can check that if 
more exact expression tp = (M^ + M^ 2 - J is used, then correction terms contribute 
already to the next-next-to-leading order. The next-to-leading order contribution 



A 2 



(34) 



can be written as the sum of partial contributions involving the propagator between 
final states, between initial and final states and between initial states, respectively: 



W (2) _ W (2) , W (2) , w (2) 



(35) 



Calculating ([34]) one obtains the following expressions for the partial contributions 

(36) 



-W {2) 
X w (f-f) 



T 



2^1- 7X^1 - 7 X 



Rh 



T X 'j 2 X X 

+ lRl T^ Rbb T^ Ra + Y Ra T^ Rbb T^x Ra + h - c - 



-w {2) 

\ (*-/) 



7^4 



rJi^j, : _ . Jin 



'i- 7 x' ■ > i- 1 x"" 

X X 

+ l 2 RbZ —Rab~ —Ra + h.C, 



(37) 



T ■ 2 »* XT ~ 1 D 



abZ ~Rn 



1 - 7 X 1 - 7 X 



'i-jX^i-jX' 

XT- 1 t XT' 1 
2 a l- 7 X aa l- 7 X 



(3* 



+ 7^a 



XT 



-1 



X 



7 



X 



l- 7 X" aa l- 7 X jRfe + 2 jR6 l- 7 X"" tl l-7X 



:Rn 



X 



R b + h.c. 



The next step is to evaluate the integrals in (|33|). This can be done by the saddle 
point method. Again, one can check that for the calculation of the cross section to the 
leading and next-to-leading orders it is enough to use the saddle point values of the 
parameters determined by the leading-order saddle point equations. These equations 
are obtained by differentiation of the expression 

- S inst (p) - S mst (p') + jW w (x ,p,p\^6), (39) 

where is given by Eq. ( j3~ID with respect to xq>, p, p', £ and 9. Physically relevant 
saddle points have x = 0, £ = and p = p', while (xo)o, £o and 6 are purely 
imaginary. Therefore, it is convenient to introduce the following notations: 

(x ) = iT, £ = ix, = -iln7. (40) 

We will denote the saddle point values of the parameters rm, x m iP an d 7 as f, x, 
p and 7 respectively. One can easily see that they are functions of e = E/E sp h and 
v = N/N sph . 

In the subsequent sections we will derive the saddle point equations explicitly, 
find their solutions and calculate the leading and next-to-leading corrections to the 
cross section in the scalar model. Here we would like to study some general features 
of these corrections. 

Suppose that the asymptotic behaviour of i? a (k) and i2&(k) for large |k| is of the 
form 

RaM ~ |k| 5 

with 5 > —1/2. This assumption is verified in concrete examples. For example, in 
the scalar theory 5 = —1/2. As it has been already explained in the Introduction, 



we will be particularly interested in the limit v — > 0. It can be shown (see Ref. fT4|| ) 
that in this limit 

X - t ~ v and 7 ~ z/ 2<5+4 . (41) 

Estimating the integrals in expression (^), one can see that the terms that do not 
vanish in the limit v — > can be written as 

W = lx<n - V In 7 - 2XS mst (p) + A [jR^T" 1 ^ + R b TR b ] , (42) 

where we found convenient to introduce the variables 

. A 



E, 9 = NX. (43) 



They are related to e and v by obvious relations 



e = -e, v = -v (44) 

K K 
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(see Eqs. 



The saddle point equations are 



dW 
dp 
dW 
d^j 
dW 

dx 
dW 







d r 



-2X^—S inst (p) + A 
op op 

-- + XRaXT-'Ra = 0, 

7 

e- A7 J R a |k|XT~ 1 J R a = 



1 R a XT~ l R a + R b TR b 



A 



jRalklXT^R. - R b w k TR b 



0. 



(45) 
(46) 
(47) 
(48) 



Here we used the definitions (j32|) , ([401). Note that the products in Eqs. ([47D and 
(|48| ) contain insertions |k| and w^. They are understood as additional factors in the 
integrand. For example, the last term in fl48|) stands for 



i dkR b (\t)w k e- w * T R b (k). 
In the limit (% — f) ~ z> — > the main contribution in the integral 
^RaXT^Ra = yj dkR a (k)e^^R a (k) 



and up to negligible corrections can be substi- 

it is easy 



and 



comes from large |k| ~ (% 

tuted by |k|, i.e. the mass can be neglected. Analyzing Eqs. 
to see that properties ( f4"TD indeed take place. 

In the next-to-leading order the terms which give non-vanishing contributions in 
the limit v — > are the following: 

= \ \^f 2 RaXT' 1 RaaXT" 1 R a + 2 1 R a XT- 1 R ab TR b + R b TR bb TR b ] +0{u). (49) 

Evaluating expressions for W and W {2 \ Eqs. (H) and ©, at the saddle point solu- 
tions of Eqs. fl4"5]) - (|48|) we obtain the functions F^'(e, u) and F^(e, u), respectively, 
in the limit v — > 0. These are precisely the leading and next-to- leading (propagator) 
corrections to the function F(e, v). 

The integrals in Eq. (]49| ) can be represented in a general form as 

i? # iO lJ R #1#2 2j R #2 = ( dkdqR #1 (k)0 1 (k)R #1#2 (k,q,e)0 2 (q)R #2 (q), (50) 



where #1, #2 = a,b (so that #1#2 = aa or a&, or 66) and the variable 9 is the angle 
between k and q. This integral transforms to 



128vr 4 / |k| 2 ci|k| / Iqfdlqli^xflkDOiflkl) 
Jo JO 



(51) 



16vr 2 



dkdq% 1 (|k|)0 1 (|k|)% 1 (|k|)5 #1#2 (|k|,|q|)% 2 (|q|)0 2 (|q|)it: #2 (|q|), 
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where the function iS#(|k|,g) was defined as 



<%(|k|,|q|) 



sin 9d6 



g#(MQ) 

16vr 2 



(52) 



for ^ = aa, ab, bb. 

Singular contributions to the propagator correction at v — > come from the first 
two terms in the r.h.s. of Eq. (pEU|). Indeed, as one can see from a simple analysis, 
due to the presence of the XT -1 factors either one of the momenta of integration or 
both are ~ m(x — t) -1 ~ m/v and are large. Hence, the asymptotic form of <S#{k, q) 
can be used. 

The leading large momentum asymptotics of the Fourier transform of the instanton 
propagator was calculated in Ref. |18[ and is equal to 



G{k, q; p) = n'(p)^(/c)^(g) In — + . . . 



(53) 



the factor n'(p) being independent of the momenta. ip^{k) is the Fourier transform 
of the zero translational mode. In a scalar theory this function can be written as 
$ M (Jfe) = k^h(k 2 )/(k 2 + m 2 ). Then 



P{s#{k, q); p) = n(p)s#(k, q) In 



s#(k, q) 



+ 



with n(p) = n'(p)h 2 (—m 2 ). From this it follows that when at least one of arguments 
is large the leading asymptotics of <S#(|k|, |q|) can be written as a sum of factorized 
terms. We obtain: 
when k, q — > oo 



5 # (|k|,|q|) = ^44/ # |k||q| 



, k i q 

In — + In — + . . . 

m m 



(54) 



when q — >• oo and k is finite 
5#(|k|,|q|) 



n 



167T 



;4/ # w k |q| 



In 



|q| 

m 



+ 



Here l aa = l bb = -1 and l ab = 1. 

Now let us analyze the first term in the r.h.s. of Eq. 
the contribution of initial states. It gives 



F^4e, v) = \Zj 2 R a XT~ l RaaXT" 1 R a 



{i-i) 
64tt 2 

A 



(55) 

which corresponds to 
(56) 



7 2 I RaXT 



-i|k| 



In 



ni 



|k| 

m 



R a XT 



-i]q| 

m 



Raj +••• 

' q 



where asymptotical formula ([54] ) was used. 
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Using Eq. (j55|) one can show that the second term in the r.h.s. of Eq. (49), 
which is due to the initial-final states, gives the following partial contribution to the 
propagator correction: 



2X^R a XT~ 1 R ah TR h 



(57) 



1 yp) \ m \m J 



R h T—Ri 

m 



+ .... 



The subscript k reminds that all factors of the corresponding integrand are taken at 
the momentum k and are integrated over k. Singularities at v — > are produced by 
the logarithmic terms: 

, |k| 1 1 
in — ~ — ~ In — . 

m X ~ T v 

Adding terms fl55D and (0) together we obtain that 



p (2) p (2) 



1 yp) y m \m ) 



x 



^RaXT-^Ra) — (RbT—Rb) 
v m J V m / q 



+ 



(5* 



In Eqs. (|5BD, (|57D and ( |55D the dots stand for terms which are finite in the limit 
z/ — ► 0. Due to the saddle equation (j48[) the expression in the square brakets in Eq. 
(|58| ) is zero. 

We would like to emphasize that our proof of the cancellation of singularities 
ln(l/V) in the propagator correction is rather general. For the proof we, essentially, 
used the general structure of the expressions for W and W^ 2 \ the leading order saddle 
point equations and the factorization property of and F^ 2 \ The latter follows 
from general formula (E3I) for the leading asymptotics of the instanton propagator. 



3 The model and the instanton propagator 



We consider the model of one component real scalar field, defined by the Minkowskian 
action 



S 



(59) 



where A > 0. The potential of the model is unbounded from below, hence the 
minimum = is metastable. Underbarrier tunelling of the initial state from this 
vacuum to the instability region and its return to the trivial vacuum is the transition 
which gives rise to the shadow process we are going to study here. 



13 



Let us consider first the case m — 0. There is a well known instanton solution in 
the massless theory given by the formula f20| , |2~l| 

<t>inst&Xo,p) = ^j=-~, — ( 60 ) 

VA (x - x ) 2 + p 2 

The solution is characterized by five parameters: four coordinates of the center of the 
instanton xq^ and its size p. Due to conformal invariance of the massless theory the 
action of the instanton does not depend on its size, 

9 (0) = <J(6 )- 167r ' 

°inst — °\H J vnst) " 



A 

In the case m^fl the mass term breaks the conformal invariance. Using standard 
scaling arguments it can be shown that there are no regular solutions of the Euclidean 
equations of motion with finite action. The decay of the vacuum <p = is dominated 
by a constrained instanton, a configuration which can be regarded as an approximate 
solution of the equations of motion. It minimizes the action under the constraint that 
the size of the configuration is p. A formalism for construction of such configurations 
and evaluation of the functional integral was developed in [p2 |. 



When pm <C 1 the constrained instanton configuration (f) c .i.( x ) behaves like the 
instanton solution (^) of the massless theory at x <C p and as a solution of the free 
massive theory, 

2^/&k pm e~ m ' x ' 



x 



vA yj\x\m \ x \ 
for x > wT 1 . The action of such configuration is 



16-7T 2 247T 2 

Si n st{p) = — —(p™) 2 



ln^ + 2^ + 1 



(61) 



4 

+ O ((pm) 4 , (pm) 4 ln(pm) 2 , (pm) 4 ln 2 (pm) 2 ) , (62) 

where Ce = 0.577216 ... is the Euler constant. For the class of constraints mentioned 
above the terms given in (^) do not depend on the explicit form of the constraint, 
whereas the correction 0((pm) 4 ) does. In our analysis we limit ourselves to the 
constraint independent order of the approximation. Therefore, all contributions to the 
propagator correction which are of the order (pm) 4 times, may be, some logarithmic 
factors or smaller are beyond the accuracy of our approximation. 

For m ^ the potential barrier separating the trivial vacuum = from the 
instability region is finite. Its height is characterized by a sphaleron solution, a static 
SO (3)-symmetric configuration satisfying the equation of motion. In Ref. it was 
found that the energy of the sphaleron is 

j j / 

E sph = K -, k = 113.4 (63) 
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and the characteristic number of particles contained in the sphaleron is given by 



N sph =j, k' = Q3 (64) 

(compare to Eqs. @ and @). 

Now let us consider the propagator in the instanton background in the massless 
theory. It is defined by the operator D x of quadratic fluctuations in the expansion of 
the action around the instanton solution (see Eq. (|])). For the massless theory this 
operator is equal to 

D x = -— + -(j)l st (x;0,p)- 



2\2 ' 



dx 2 ^ 2 ms ' dx 2 (p 2 + x 

It can be easily seen that it possesses five zero modes i[)a(x) (A = 1,2,3,4,5) cor- 
responding to the translational invariance and the scale invariance of the massless 
theory. The zero modes can be obtained by differentiation of the instanton solution 
with respect to the parameters Xq and p: 

d 

4>A ~ -W7~(pinst(x; X , p) \x =0, Cm = («o)m> Cb = P- 

In accordance with the general theory the propagator in the instanton background is 
the inverse of D x on the subspace of functions orthogonal to functions /U(x) which 
satisfy the only condition that the matrix 



^ab = J dxip A (x)f B (x) 
is invertible [p3|| . Let us denote such propagator by Gf(x,y). It satisfies the equation 
D x G f {x, y) =6(x-y)-Y, J : a{x)Q^ b(x) (65) 

A 

and the orthogonality constraints 

dxf A (x)G f (x,y) = = / dyG f (x,y)f B (y). (66) 



The lower index / reminds that the instanton propagator satisfies the constraints 
defined by f(x). The r.h.s. of Eq. ( pB] ) is the projector onto the subspace orthogonal 
to the functions Ja{x)- Thus, there is an ambiguity in the definition of the propagator 
due to the choice of the functions f A (x), but physical results, of course, do not depend 
on it. In the next section, following the ideas of Refs. |12|, |2^| , we will use the freedom 



of choosing constraints (p6|) to eliminate the leading asymptotics of the instanton 
propagator and simplify the analysis of the initial-state corrections. 
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Let us consider first a particularly simple and natural choice of the functions 
f A (x), namely 

f A {x) = w(x)ip A (x), 

where the weight function 

w(x) 



(p 2 + X 2 ) 2 

Zero modes i/)a(x), orthonormal with respect to this weight function, are equal to 
, ( , P /l5 Ix^p p /l5 x 2 - p 2 

In this case VLab = &ab- Allowing some abuse of notation we denote the propagator 
orthogonal to the zero modes themselves by G^(x,y). It satisfies the equation 

/ d 2 24p 2 \ 5 

+ < 2 i 2^2 G t( X ' V) = S ( X ~ V) ~ W ( X ) 12 ^A(x)lpA{x). (68) 



dx 2 (p 2 + x 



A=l 



Making use of the fact that the symmetry group of the (massless) theory is the 
large conformal group SO (5) one can reformulate the theory as a free theory on 
the four-dimensional sphere S 4 ]TD|, ^TJ. The relation between points on R 4 and 
corresponding points on S A is given by the standard stereographic projection. The 
instanton propagator G^(x,x') on R A is related to the (free) propagator Ify (£,£') on 



S by the formula |]T9 



where the points £ and £' on the sphere correspond to the points x and x' in the 
Euclidean space respectively In fact, it can be shown that the propagator depends 
only on the geodesic distance s(£, £') between the points. It is convenient to introduce 
the function 

d(£,0 = ^(l-coss(£,0)- 
As a function of the coordinates of the Euclidean space it is equal to 

The propagator K^(d) satisfies the equation 

(-a c - 4) j^(d(e,o) = - E movUo, (71) 
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where Ag is the Laplacian, is the <5-function on S 4 , and the orthor normalized 

zero-modes ^a(0 of the operator in the l.h.s. of Eq. (|Tl) ) are related to the zero 
modes i/)a(x), Eq. (p7j), by 

If the sphere S* 4 is considered as being embedded into the 5-dimensional Euclidean 



space with coordinates {za} (A = 1,2,3,4,5) then ipA = J15/87t 2 Za- 

There are a few ways to find the instanton propagator K^(d) satisfying Eq. (|7T|). 
Perhaps the simplest one is to start with the expansion 

oo I r 

Z=0{m} A l 4 

The spherical harmonics Yi, m (£) on •S' 4 are eigenfunctions of the Laplace operator with 
the eigenvalues A;: 

-A^, m (0 = A,Y,, m (0, 

A, = /(Z + 3), / = 0,1,2,..., 

where Z is the total momentum and m denotes the set of orbital momentum numbers. 
Using summation formulas the explicit expression for the instanton propagator on S* 4 



was obtained IS 



111 2 



1 43 56 ' 
— - — 3 In d h 6d In d -\ d 

2d 5 5 



Using Eq. (|69|) we arrive at the final expression for the instanton propagator: 

G^(x,y) - 



2tt 2 (p 2 + x 2 )(p 2 + y 2 ) [2d(x,y) 

43 56 1 

- 3\nd(x,y) - — + 6d(x,y)\nd(x,y) + —d(x,y)> , (72) 

where the function d(x, y) is given by Eq. (|7Dp. 

When points x and x', or the corresponding points £ and approach each other 
the propagator G^(x,x') is singular, the leading singularity being given by the first 
term in the curly brackets in the l.h.s. of Eq. ([73). This is precisely the singularity 
of the free propagator of the massless scalar theory in the four-dimensional Euclidean 
space. Indeed, this propagator, denoted here by G (x — x'), satisfies the equation 

d 2 

G (x — x') = 8(x — x') 



dxl 
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and is equal to 

Go{x ~ ' x ' ] ,= -wj^hf (73) 

The corresponding propagator K (£, £') on S 4 , related to G (x — x') by the formula 
similar to (|69"D , satisfies the following equation on the four-dimensional sphere: 

(-Af + 2) Jfotf, O = O 
(compare with Eq. (|71|)) and is equal to 

/ 1 1 
M&O = 8^2 2 d(f,f')' 

This gives the first term in the curly brackets in Eq. ( ff^) . 

Let us discuss briefly the structure of the expression for the instanton propagator. 
A detailed analysis shows [19] that G^(x,x') can be written as the following infinite 
sum: 

oo 

G^(x,x') = G (x — x') + G n (x, x ) 

n=l 

A r . „ . , l2 



= G (x - x') + - J dyG (x - y)(f) inst (y)G Q (y - x ) + . . . (74) 

This series has an interpretation in terms of Feynman diagrams. The term G n (x,x') 
is represented as the diagram consisting of a line, corresponding to the free scalar 



propagator (ff3|), with n insertions (j)f nst - The first term in Eq. (|74]) is just the free 
propagator without insertions. As we have already discussed, it gives the leading 
singularity of the full propagator G^(x,x') when x — > x'. The integral term, written 
down in Eq. (|74]) , is the Gi(x,x') term. Calculating it we get 



2tt 2 (p 2 + x 2 )(p 2 + x' 2 ) — p 2 (x — x') 2 (p 2 + x 2 )(p 2 + X 2 ) 

This term has the logarithmic singularity when x — ► a;'. This is precisely the sublead- 
ing logarithmic singularity of the exact expression (|72|). 

We finish this section by presenting a relation between the propagator Gf(x,y), 
satisfying a general constraint (|BUD, and the propagator G^(x,y). They are related 
as follows: 

Gf{x,y) = G^(x,y) - (J dzG,p(x,z)f A {z)^VL A l B ip B {y) 
- iI)a(x)(Pab) J dz fB(z)G i> (z,y) 

+ ^(^(^^(/^^(^^(z^O/c^'))^^^)- (75) 

In the next section we will see that with the help of this formula and by an appro- 
priate choice of the function f(x) one can modify the asymptotics of the instanton 
propagator. 
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The high energy asymptotics of the instanton 
propagator 



The Fourier transform of the instanton propagator is defined by Eq. (|22|) . In princi- 
ple, using the exact expression (|72|) the function G^(p,q) can be obtained by direct 
calculation. We did not find the complete analytical expression for it, instead we 
derived the asymptotic formula for the Fourier transform of the instanton propagator 
in the regime when p 2 , q 2 are fixed and s = (p + q) 2 — > oo. The growing terms of the 
asymptotics are given by 



G^(p,q) 



where 



16tt 2 
p 2 q 2 



sp 2 ln( S p 2 )iTi(p, q ) + (sp 2 )U 2 (p, q) + \n(sp 2 )U 3 (p, q) + . . . , (76) 



niO,g) 

n 2 (p,g) 
n 3 (p,g) 



S 1 {j>p)S 1 {qp) 

1 



+ 



51(pp) 

9 

-S 2 (pp) - 



15 
9 



ln2)5 1 (pp)5i(gp) 
27 3 



2 2 



27 



' 2 2 

rP P 



Si(pp) 



Si(gp) 
^i(gp) 



-S 2 (pp)S 2 (qp) 



Here <S n (z) is defined by S n (z) = z n K n (z), where K n (z) is the modified Bessel func- 
tion. Using the explicit expressions for the normalized translational zero modes Eq. 
([67]) , one can easily see that the product of their Fourier transforms ip^p^^q) 
Thus, the first two terms of the asymptotic expansion (|76|) can be written as 



p 2 s. 



G^{p,q) 



1 



5p 2 



5p 2 



1 



\n(p 2 s)^(p)^(q) - — (C E - — - In 2 %{p)%(q) +... (77) 
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The leading term of the asymptotics of the propagator in the instanton background 
was calculated in Ref. |0| and is given by Eq. (|53|) . This result is in complete 
agreement with the first term in Eq. ([77]). 

In Ref. 0] Mueller proposed the idea to use the ambiguity in the choice of the 



function Ja(x) in order to cancel the two leading terms in the asymptotics of the 
propagator Gf(p,q). If this can be done, then the propagator contribution and loop 
contributions of the initial state corrections vanish. As a consequence, such correc- 
tions do not exponentiate, i.e., do not give contributions to the function F(e,v). 
Moreover, in this case the initial-final state corrections can be described semiclas- 
sically. Namely, the effect of the initial-state lines can be taken into account by 
substituting the instanton by a new field configuration which is a particular solution 
to the classical equation of motion with an external source (see Ref. [I2| for details). 
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In the rest of this section we discuss how the functions Ja(x) that provide vanishing 
of the two leading terms in Eq. ( |76|) (or Eq. ([H])) can be chosen. For this we 
essentially repeat the arguments of Ref. []I2"| . The propagator constraint for 
which the vanishing takes place turns out to be not relativistically covariant. Let p\ 
and P2 be the arguments of the Fourier transform of the propagator. We choose a 
coordinate system such that the components pij = p 2 j = for j = 2, 3, whereas 

Pi+P = P2P > 1 

and p\ and p\ are fixed. Here the ± components of the momenta are defined by 

0j)o ± (Pj)l 



V2 



Then (pi,p 2 )p 2 ~ P1+P2-P 2 1- Only the components f^x), corresponding to 
translations, play a role. The Fourier transforms f^{p) of the functions ffj,(x), defining 
the required propagator constraint, can be chosen in the following way: 

Up) = S(p + - M)5(p_ + M)f,(p 2 ,p 3 ), 

where M is an arbitrary parameter of the dimension of mass. Substituting such 



functions into Eq. ( [75]) one finds after some calculation that indeed the s In s- and 



s-terms of the asymptotics ([fg) cancel and 

Gf(pi,p 2 ) = ^y,{Pl)Qnv^ V {p2) + • • • , 

where the 4x4 constant matrix is equal to 

^ = -^f^J^f J AiA 2 /.(-gi)^(gi)ln^^(g 2 )/ T (-g 2 ) (n T )-\ 

Using the freedom of choosing the function f ^ one can make the constant real sym- 
metric matrix Q^ v equal to zero. We would like to stress that the knowledge of the 
exact expression for the instanton propagator allows us to get the explicit formula 
for the matric Q^. This is in contrast with the case of the electroweak theory, where 
only a general structure of the analogous matrix can be derived 



5 Residue of the instanton propagator 

As we have seen in Sect. |2] for the perturbative calculations of the function F(e, u) 
the on-mass-shell residue of the instanton solution is needed. Its definition is given 
by Eq. (0). In the massless case the Fourier transform of the instanton solution is 
equal to 

4> mst (p;0,p) = ^=^-K 1 (p\p\), (78) 
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where R = IQy/Zir 2 , K x (z) is the modified Bessel function, and 

P 



I(p) = -^=R . (79) 



Using definitions flIE|), ( |TTD we obtain that 



1 p 16^3 



7T 



2 



«= & < k) = (W' w ^(W' (80) 

In massive theory one should take the Fourier transform of the constrained instanton 
solution. Later we will show that, in fact, within the approximation considered here 
it is enough to take the residue (|79|) of the massless instanton. Corrections due to 
non-zero mass are of the order 0(p 4 m A ), i.e. of the order of terms already neglected 
in (pip. However, in Eqs. fl8U| ) we should use = Vk 2 + m 2 , i.e. the expression for 
the energy of the massive theory. 

To calculate the next-to-leading correction to the function F(e, v) we need the 
expression for the double on-mass-shell residue of the instanton propagator. The 
propagator residue is defined by Eq. (|2~6|). In the massles theory we have 



^^^^(M, (81) 

where # = aa, ab, bb and the function s^(k, q) is the s- variable for the corresponding 
particles on the mass shell: 

*2(k,q) = 4?(k,q) = -aSfrq) = -2(|k||q| - kq). (82) 



The functions P(p 2 s) and P(s; p) in Eq. (|2~6D are related in the following way: 
P(s; p) = p 2 P(p 2 s). However, in the calculation of the next-to-leading order cor- 
rections due to non-zero mass must be taken into account. It turns out that within 
the accuracy set by Eq. (|H]) it is enough to use functions -R#(k, q) defined by the 
relation 



where, as in the case of Eq. (pT|) , P is the residue of the instanton propagator of the 
massless theory, whereas the energy and the s- variables are taken for the massive one. 
Namely, = Vk 2 + m 2 and s#(k, q) are given by Eqs. (24) , (|25|) . The motivation 
for such procedure of calculation is discussed in Sect. ^. 

Note, however, that there is an ambiguity in representation (^). One can equally 
well substitute s^P (k, q) in (|81 ) with any other function of k and q which coincides 

with s^(k, q) for k 2 = q 2 = 0. For example, one can use the expressions coming 
from s = 2(k, q), i.e. 

s aa (k,q) = s bb (k,q) = -2(^ u ^ q - kq), (84) 
s ab (k,q) = 2(^ k u; q -kq). (85) 
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We will see later that the effect of this ambiguity is negligible for the final result. 

The exact expression for the function P(p 2 s) was obtained in Ref. |25|. Its calcu- 
lation is a rather tedious although straightforward procedure, and we do not give the 
details here. Instead we would like to discuss some general features of the computation 
before presenting the answer. 

The terms in expression (|72|) give contributions to the residue which can be divided 
in four classes. 

1) The first term in the curly brackets in Eq. (72) gives rise to the free propagator 
as it was already explained. Its contribution to the Fourier transform G^(k,q) is 
proportional to 

2(2tt) 4 



k 2 



-5(k + q). 



This describes free motion of the particle not interacting with the instanton and is 
irrelevant for our problem. 

2) There are factorizable terms of the form gi(x)g2{y), where gi(x)'s are propor- 
tional to expressions like 

x n x n ln(p 2 + x 2 ) 

(p 2 + x 2 ) 1 ° r (p 2 + x 2 ) 1 { } 

with some integer n and I. Their contributions to the momentum-space propagator 
are of the form gi{k 2 )g2{q 2 )- These are s-independent contributions to the function 
P(p 2 s). 

3) The next group of terms is of the form (xy) g\(x) g2(y) with g^ of the form (|86|). 
Calculating the momentum-space propagator we get 

, ik * +i ™(x,y) gi (x)g 2 (y) = - — — 



ikx+iqy 



d d 



Unu<t) 



dk» dq>* 
-^qyg'^yg^q 2 



gi(x)g 2 (y) 



dkv dq^ 

This gives a contribution to the residue proportional to s#(k, q). 

4) The last group consists of terms of the form (xy)ln(x — y) 2 gi{x)g2{y) and 
ln(x — y) 2 gi(x)g 2 (y). Carrying out the computations one can show that they lead 
to terms proportional to s#lns# and lns#, respectively, in the expression for the 
residue of the instanton propagator. 

Finally, the exact expression for the function P(z) is given by 



P(z) 



16?H 



-zln- + -z { C E - — 



2 4 



If) 



Below this result will be used for the calculation of the next-to-leading correction 
to the function F(e, v). It is convenient to write Eq. fl87|) as 



P(z) 
a i 



3/4, 



Q' 2 



zln- + 2 
4 

= 3/2. 



C E - 



1 

15 



— a 2 



In - + 2 [C E + 



43 
30 
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This form of the residue allows to trace the origin of various contributions to the final 
answer for the propagator correction. 

From Eq. ([H]) we see that the important ingredient of the contributions to the 



propagator correction is the function <S#(|k|, |q|), defined by Eq. fl52|) . With the 
residue P(s) given by expression (|S8|) <S#(|k|, |q|) can be easily calculated. As it was 
already discussed in Sect. |2|, when one of the arguments is large the asymptotics of 
this function can be represented as a sum of factorized terms. In the model under 
consideration we obtain 
when k, q 



OG 



<S#(|k|, |q|) ~ 4/ # |k||q|p 2 a! 



ln(|k|p) + ln(|q|p) + 2 In 2 



C A 



i9) 



when q — > oo and k is finite 



«S # (|k|, |q|) ~ 4/ # w k gp 2 ai 



ln(|q|p) + ln(|k|p) + In 2 + C A + M#j2 (|k|, oo) 



(90) 

where l aa = = —1, l a (, = 1 and M #)2 (|k|,oo) is given by Eqs. ( |A2|) in the 

Appendix. Let us recall that the leading terms in formulas (|S9"D, ( PD| ) were obtained in 
Sect. H for a general model from the leading asymptotics of the instanton propagator 
ofRef. ill 



6 Leading and propagator corrections: case v —> 

In this section, firstly, we study the saddle point solutions that give dominant contri- 
bution to integral ( pB]) in the limit v — > 0. Secondly, with this solution we calculate 
the leading and next-to- leading order contributions to F(e,v). Finally, we give an 
illustration of cancellation of terms singular in v in the case of the concrete model 

For this model the leading contribution to F(e, v) comes from the term 

W = A8n 2 p 2 m 2 \n(C p 2 m 2 ) + ixm - ulwy (91) 
+ 192vp 2 [J„(7, r, x) + iMl, X~T,x) + 2 7^o(7, X, x)] 



(see Eqs. (p0|), (pT|) and ( pi]) ). Here we introduced the notations 

lnC= -In4 + 2C £ + 1 « 0.768137. .. , v = 192n 2 

and 

Join^x) = t~ / — 

47r J Wk 1 — 7e 



-«>kX 



Recall that e and z/ are defined by Eqs. ( pE3"D and are related to e = E/E sph and 
z/ = N/N sph by relations ( f44|) with k and given by Eqs. (|63|), (|64|). 
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In accordance with estimates ([11]) we expect that in our model 

X — f ~ v, and 7 ~ zA 



(92) 



It will be shown shortly that this is ineed the case. In the leading approximation in 
v integrals J in Eq. ( pi]) can be taken at 7 = 0. Then they can be easily calculated, 
and one gets 

$(rm) = -Lj(Q,r,x) = ^iH. (93) 



TTfl 



Taking into account (92) we find that in the leading order in v the following expression 
for W can be used: 



W 



48ir p m ln(Cp m ) + exm — DYn."{ 



(94) 



+ vp 2 m 2 



$(rm) + 



7 



+ 0{v 3 



Then the saddle point equations 




m 2 (x — t) 2 
8|) take the form 



■W 

dp 2 

d_ 

dx 
d_ 



48vr 2 m 2 In (p 2 m 2 Ce) + 
1 



vm 



$(rm) 



7 



m?(x ~ r ) 2 



z> 



7 (X - r ) 



me — 2-up 



7 



■up 



(x-r) 3 ' 



m 3 $ mr + 



(95) 
(96) 
(97) 
(98) 



(x-r) 3 . 

as in Sect. || let us denote the saddle point solutions for the dimensionless parameters 
mr, mx, Tap and 7 as f, x, p and 7 respectively. We are looking for these saddle 
point solutions as functions of e and v. Here and below allowing an abuse of notation 
we will use the same letters for functions of e, v and e, v. This will not lead to a 
confusion. 
From Eqs. 



- fl9"8" ) in the leading order at v — > we get 
1 I 



P 2 



7 
X 



V 

f + 2- 

e 



The function r(e, z/) is a solution of the equation 



In 



eCe 



+ 4 



7/ 



(99) 

(100) 
(101) 

(102) 
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Figure 1: Plot of function e(ro) for v — 0. 

This equation is, of course, a consequence of the system of the saddle point equations 
([95]) - (p8|). Its solution can be obtained as an expansion in v. We write 



r.2\ 



f(e, z/) = f (e) + £f 1 (e) + 0(z? 

After substituting this expression into Eq. (|102|) one finds that the coefficient f\ is 
equal to 

- _ 4 ($'(f )) 
1 e4($'(f ))-$"(f )' 
whereas the function fo(e) is determined by the equation 

ln bH)) +4 * (foHO - < 103 » 

Note that the solutions (|100D , ( gQ| ) are in agreement with the behaviour (^) assumed 
above. The behaviour of the solution 7o(e) can be analyzed by writing first equation 



(|T03|) as 

studying the function e(ro) (or e(r ) = e(r ) / k) and then inverting it. The plots of 
functions e(r ) and f (e) are given in Fig. [I] and Fig. |2|. 

Let us write the exponent F(e, z/) of the cross section (d) as 

F(e, v) = -32tt 2 + F«(e, v) + F^ 2 \e, v) + . . . , 
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Figure 2: Plot of function fo(e) for v = 0. 



where F^(e,u) and F( 2 )(e, z/) are given by W and W^ 2 \ respectively, evaluated at 
the saddle point solution. Substituting solutions (|9"9"| ) - ( |101| ) into Eq. ( |9~3] ) we obtain 
that 

" ~ / \ 1 
r (e) + 



F (D, 



e, i/i 



4$'(f (e)) 



(104) 



For < f < 1.28 the function e(fo) grows monotonically from e = to e « 0.59 
(see Fig. [1]). Hence, for the interval < e < 0.59 there exists the inverse function fo(e) 
with the property t (0) = 0. This is also the interval of energies for which the saddle 
point solution at v = can be found. The plot in Fig. ^] shows that for < e < 1.59 
the values of T (e) in general are not small. Since f is the saddle point solution for 
rm, this observation suggests that the mass corrections cannot be neglected. For 
< e < 0.1 we enter the regime when f < 1 and the expressions for the saddle point 
solution and Eq. ( |104j ) simplify considerably. In particular, $(f ) ~ 1/tq for f <^ 1. 
In this regime corrections due to non-zero mass can be neglected. 

At very small energies, namely when 



mini 



< 1, 



(105) 



we enter the regime when analytical expressions for the saddle point solutions can be 
obtained. The function fo(e) is small in this limit. In the leading approximation Eq. 
(11031) takes the form 



4 

In e + -ttj + 5 In r + In 



7n 



1 + ... = 0, 
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where the dots stand for terms vanishing as e — > 0. It is easy to obtain the asymptotic 
form of the solution of this equation in for the variable l/f (e). We get 

j 115 1 

- In — | — In In — h In v — 4Cr — 1 + decreasing terms. 
4 e 8 e & 



The first two leading terms of the function T (e) are given by 

, x 2 5 lnlni 



H 2 (m|) 3/2 



+ . . . 



(106) 



Then in the leading order in energy solutions (^) - (|101| ) become 



P 2 



48tt 2 ( ln l) 3/2 ' 



7 



zA 3 / 1\ 3 / 2 



(107) 



In this regime the function F^(e, u) is equal to 



F (D( 



ln± 



1 + 



' In In - 

~hTr 



(108) 



This is result fl6|) quoted in the Introduction. 

The propagator correction F^ 2 \e, u) in the limit i/ — » is obtained by evaluating 

) - (|101| ). In Sect. we showed 



contributions (^) - (p8|) at the saddle point solution 
in a rather general context that singular terms ln(l/i/) appear in the contributions 
involving initial particles and proved that they cancel each other in the final result. 
Let us discuss now the exact expression for the propagator correction for the (— A0 4 )- 
model in the limit v — > 0. The result for the partial contributions is given by Eqs. 
(|A3|) - (|A8|) in the Appendix. We see that singular terms ln(l/z>) appear only in 

(2) (2) 

the functions F^_^ (e, v) and of Fh_^\ (e, v) which involve initial states and which are 
proportional to a,\. Hence, according to Eq. (p8|), the singular terms are due to slns- 
and s-terms in the instanton propagator. It is clear now that, if following the method 



of Ref. 



reviewed in Sect. 



one chooses the propagator constraint in such way 
that these terms are absent in the residue, then the contributions of the initial states 
are zero in the limit v — > 0. Correspondingly, in this case no singular terms appear. 



As it follows from these formulas and the saddle point solution ( |106|) , in the limit 

") behave as 



of small energies the corrections F^* 



and of F^l f) ( 



H) 



3/2 ' 



(2) 

The contribution i^_«(e, v) of the final states contains a part ~ a± and a part ~ a^- 
The former is due to sins- and s-terms in the instanton propagator residue, the 



27 



latter is due to Ins-terms and s-independent terms of it (see Eq. (|88|). In the limit 
of small energies F^_^(e, u) behaves as e 2 , thus giving the dominant contribution to 
the propagator correction. 

As an illustration let us demonstrate the cancellation of the singular terms explic- 



itly in our model. Summing up corrections (|A4j) and (A6) we get 
F$l l) M + F$ ) _ l) M = ?>2a 1 v~p« 



X 



7 



2(x " ^ ) 3 

7 
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(x - fy 
p 2 



71 , P 2 
30 + n (x-r) 2 



$'(f) ln^-r 



7 



r 2 2(x-f) 3 



$'(f)(21n2-^) + 



7 



(X 



The singular term is ln(p 2 /(x — f) 2 ) ~ ln(l/P) in the second line. The coefficient in 
front of it vanishes exactly due to the saddle point equation fl98l) . 

Finally, let us present the result for the complete propagator correction to the ex- 
ponent of the multiparticle cross section in the limit v — > 0. Summing up expressions 
(A4), (A6) and (|A3) for the partial corrections, we obtain 



F^(e, 0) 



a 2 ~ 2 4$ 2 (f) 
~ e ($'(f)) 2 



8«i e 



43 
15 



f 2 



?3 



■u 2 $'(r 
($'(f)) 2 



- - In 2 



ln2-2C E -4$(f) - In— + 



4 ' $ 2 (f) 



($'(f)) 2 



(109) 



ln2-^-2C £ -4$(f)-21n^ 



The expression for the propagator correction simplifies in the regime f < 1. 
In this case, using Eqs.( |AlO| ), ( |A1 1| ) and ( |A12| ) or calculating the limit of small fo 
directly from Eq. (|109|) , one gets 



F^(e,0) = ^e 2 



v 



f 2 



l + ^(81n^ + 8C s -^)+O(f 4 ) 



43 \ 



15/ 



(110) 



4ai 



e To) 3 -0(f 2 ,f 2 lnf ). 



The first line is the leading term at small energies: 

F( 2 )(e,0) = ^e 2 - 1 



2 2 



'Ill) 



v 32tt 2 

As it was discussed above, this is the contribution from the propagator between final- 
final state. The first correction to expression (|111|) is given by the f 2 lnf term in Eq. 
( |110|) . In the regime of very small energies we have 

2i mini 
r In r 



In 
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Figure 3: Region of values of e and v for which the saddle point solution exists. 



By comparing the propagator correction ( 111 ) with the leading order correction ( |108| ) 
we see that the actual expansion parameter at small energies is indeed eJhx{l/e), as 
it was argued in Ref. fI3|. 



7 Propagator correction in general case 

For arbitrary e and v the saddle point equations are obtained by calculating of the 
derivatives of the general expression W, Eq. (0), with respect to p 2 , 7, x an d t. In 
this case we found the saddle point solutions p 2 (e, u), f(e, u), x(e, v), j(e, v) for imp) 2 , 
mr, mx, 7, respectively, and computed the functions fw(e, u), F^(e, v) numerically. 
It turned out that the saddle point solution exists only for the region in the (e, v)- 
plane which is presented in Fig. [| It lies inside the rectangle < e < e max = 0.59 
and < v < v max = 0.25. For points very close to the axis v — our numerical 
computations fail. Extrapolating these numerical results to v = we found good 
agreement with the analytical results for v = obtained in Sect. |6|. 

We performed the numerical analysis for the whole region, plotted in Fig. |3|. For 
the presentation of the results it is convenient to introduce the function jF(e, z/), 
related to F(e, v) by 

F(e,u) = -327r 2 .F(e,z/), (112) 
and its leading and next-to-leading approximations: 

FU(e,u) FW(e,v) + FM(e,v) 
*k(e, v)=l- and T 2 {e, u) = 1 ^ . 



29 



All these functions are normalized by the conditions .F(0, v) = .Fi(0, v) = -7-2 (0, v) = 
1. 

Before presenting our results we would like to mention that the complete function 
F(e, v) in the range 0.4 < e < 3.5 and 0.25 < v < 1 was computed in Ref. f9j. The 
computation was performed by solving a certain classical boundary value problem 
on the lattice. With the size of the lattice used in the numerical calculation in Ref. 
0, the authors did not obtain data for smaller e and v except for the line of points 
corresponding to the periodic instanton solutions PB[| . This line is directed from the 
zero energy instanton (e = v = 0) to the sphaleron (e = v = 1). 

The left upper bound of the region in Fig. ^| is precisely the line v = v p (e) of 
values of e and v of periodic instantons. For them r(e, v) = x(e, z/)/2. The leading 
order and the next-to-leading order approximations for jF(e, v) at the line (e, u p (e)) are 
shown in Fig. £|. One can see that in the whole range of calculation these curves are 
close to each other. They can be compared with the complete function jF(e, f p (e)) for 
the line of periodic instantons, computed numerically in Ref. @ and also plotted in 
Fig. The comparison shows that our perturbative results do not differ significantly 
from the exact ones for e < 0.25 and v < 0.2. These values can be regarded as 
a rough estimate of the range of validity of the perturbative calculations up to the 
next-to-leading order. 

Lines of constant T%{e, v) are plotted in Fig. || Within the range of (e, v) in Fig. || 
the minimal value of T{e 1 u) is a bit below T = 0.77. This gives still a considerable 
suppression of the multiparticle cross section of instanton induced processes. 

As it was explained above, our main interest is to study the cross section for 
shadow processes with a few initial particles. According to conjecture fllPp the points, 
where the lines of constant JF 2 cross the v — axis, are of particular interest. Ex- 
trapolating our numerical results to v = 0, one obtains, for example, JF 2 (e, 0) = 0.95 
at e = 0.180, JF 2 (e, 0) = 0.85 at e = 0.492. These values, of course, can be obtained 
by using the exact formulas for v — from Sect. |6|. We would like to mention that 
in the region of (e,v) which we studied the next-to-leading approximation jF 2 (e, z/) 
is quite small comparing to the leading order J-\{e, v). The difference between these 
functions does not exceed 10~ 2 . The curves in Fig. § end at the line formed by saddle 
points of the periodic instanton solutions. 

In Sect. [| we mentioned that various expressions for the functions s#(k, q) can be 
used in the definition of the residue of the instanton propagator, Eq. (B3f). Examples 



of such functions are given by Eqs. (P4|) , (f25|) and (]84|) , (|85|) . We performed the cal- 
culation of the leading and next-to-leading order approximations for these particular 
choices of functions and for a wide range of values of e and v. We found that the 
difference is very tiny. Thus, 

^ 2 (0.2,i/p(0.2))| Wrt -^(0.2,^(0.2))!^ w 5 ■ lO" 4 . 
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epsilon 



Figure 4: Plots of J-\{e, v p {e)) (the lower curve) and Ti{e y v v {e)) as functions of e for 
the periodic instanton. For comparison the complete function J-(e, v p {c)) (the longer 
curve), calculated in Ref. ||, is also plotted. 
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Figure 5: Lines of constant JF 2 (e, v) in the (e, z/)-plane. Numbers near the lines 
indicate the value of JF 2 , "p" labels the line of periodic instanton solutions. 
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8 Accuracy of the approximation 



As it was explained in Sect. || in principle the constrained instanton configuration and 
the action of this configuration depend on the constraints imposed. In Eq. (51) we 
calculated only the leading, constraint independent part of the action. This is main 
restriction to the accuracy which can be achieved in our calculations. Corrections 
due to the form of the constraint after evaluation at the saddle point solution are of 
the order p 4 , p 2 lnp or p 2 In 2 p (see Eq. (|6l|) ). It turns out that within this accuracy 
main terms of the leading and propagator contributions to F(e, v) can be calculated. 
For this one has to follow a certain scheme of approximations, which we discuss in 
the present section. 

Let us first obtain some simple estimates of the leading and next-to-leading cor- 
rections at v — > 0. For the case when mr <C 1 the dominating term of the leading 
correction is given by p 2 Jo(7, t, x) ~ P 2 / r2 - From the analysis of the next-to-leading 
one can show that the main term is F^_^, which is due to the instanton prop- 
agators between the final states. For small mr the dominant contribution to the 
integrals in Eq. (|36|) are given by the soft momenta |k| ~ |q| ~ 1/r. For such 
momenta the ln(p 2 s)-term in the expression (|5^) for the residue of the instanton 
propagator plays the main role. Putting all the factors together one can see that 
_p( 2 ) ^ F(j_^ ~ (p/r) 4 ln(p/r) 2 . To summarize, in this regime 

^ (2) M)~4 in 4- ( n3 ) 



These estimates are in agreement with results ( |104j ), (|110|) and (|A7|) . To see this one 



er ~ p 2 r$'(r) ~ ^ 



should use formula (^) for p 2 and take into account that for r <^ 1 

I 

f 2 

in Eq. (gDg), and 

n 2 n 4 n 2 

r A r 4 T A 

in Eq. (|A7|). For very small energies from solutions ( |106| ), ( |107| ) one gets 

t 2 ./ m i r 4 r 2 



in accordance with results ( |108| ) and (|111 



To have numerical estimates of the range of values of these characteristic terms we 
studied the saddle point solution p 2 (e,u) and the ratio p 2 (e, z/)/f 2 (e, v) numerically 
for a wide range of e and v. In Fig. |6| and Fig. [7] we present the plots of these functions 
for the most interesting case v = 0. 
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Figure 6: Plot of the function p 2 (e, 0). 




For the procedure of calculation to be consistent one has to check, first of all, 
that p 4 , p 4 lnp and p 4 ln 2 p are smaller than the propagator correction F^ 2 '(e, v). We 
checked numerically that this is indeed the case for a wide range of e and v. In general, 
since, as it was mentioned in Sect. [7|, \T\ — T?\ < 1CT 3 , the propagator correction 
< 0.3, whereas p 4 < 10" 4 . Thus, at the point e = 0.2, v p (0.2) = 0.17 on the line 
of periodic instantons F^(0.2, 0.17) = 0.28, whereas p 2 (0.2,0.17) « 0.02. Another 
example: for v = F( 2 )(0.33,0) « 0.3, whereas p 2 (0.33,0) « 0.01. From Fig. §) 
and Fig. [?| one can see that the corrections p 4 , p 4 In p and p 4 In 2 p are smaller than 
F( 2 )(e, v) for the regime when r < 1 (see Eq. ( |113[ ). This conclusion is also confirmed 



for v = in the limit of very small energies. Using formulas (|106|) , ( |107|) one can see 
that 

f 2 f 2 f 2 



p ~ - — —3, p lnp ~ - — —2, p In p 



In 



Recall that the propagator correction F^(e, 0) ~ e 2 in this regime (see Eqs. ( |1 1 1| ) 
and (|114|) ). To give an illustration we presented the plot of the function p(e, 0) in 
Fig- § 

It turns out that with the accuracy set above it is enough to use the residue of 
the instanton solution I(p) and the double on-mass-shell residue of the instanton 
propagator P(k, q) of the massless theory. Let us justify this point. 

Let <fic.i.{x\p) be the constrained instanton solution of the size p in the massive 
theory. From general arguments one can see that its Fourier transform can be written 
as 

7 /, n \ g(k 2 p 2 ,m 2 p 2 ) 

<t> c ,.(k-Ap)=P k2 + m2 ■ 

Expanding the function g(k 2 p 2 ,m 2 p 2 ) in powers of m 2 we write 

g(k 2 p 2 , m 2 p 2 ) = g (k 2 p 2 ) + (pm) 2 g 2 (k 2 p 2 ) + ... 

Logarithmic terms of the form (pm) 2 ln(pm) may also appear, however we will not 
write them explicitly assuming that they are roughly of the order (pm) 2 . The residue 
Ic.i{p) = {k 2 + m 2 )(j) c .i(k] p)\k =iw k is equal to 

h,{p) = pgo(0) + p(pm) 2 (<? 2 (0) - g' (0)) + ■■■ 

Since the constrained instanton solution reproduces the instanton solution (p inst (x; 0, p), 
Eq. (^3), of the massless theory in the limit m — > 0, the Fourier transform of (pi ns t 
and its residue can be written as 



4>inst(k; 0, p) = p — , I(p) = pg {0). 



goipp 2 " 
k 2 

These formulas are in accordance with exact expressions (|78|) and ([791). 
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Typical terms of the leading order correction are of the form 

dk 



where 



l cA . (p)K(w k T, w k x)I c .i. (p) , (115) 



K(w k r, w^x) = z z— (116) 



(see Eq. (|91|)). By a simple analysis one can show that 

dk 



1 



I c .i.{p)K(W] s T, WkX)Ic.i.(p) 

w k 



f dk ( P 2 \ 

= / — hnst{p)K{w] s T,W] / oi)Iinst{p) + O mp 2 — . 

J W k \ T 2 J 

We checked that the ratio p 4 /f 2 is small comparing to F^(e, u) for a wide range of 
e and v. This can also be seen from the plots in Fig. |B|) and Fig. 0. In the regime of 
very small energies from Eqs. ( |106| ), ( |107| ) one gets 

whereas F^(e, 0) ~ e 2 . 

However, if the expression % = |k| of the massless theory is used in ( |115|) , then 



by simple estimates one can see that 
dk w s _ , , , rr , f dk 



I —I(p)K(w k r, w kX )I(p) ~ / £l(p)K{\k\T, \k\ X )I(p) + 4 • 0(m 2 r 2 ). (117) 



Since the saddle point solution f(e, v) > 1 for a range of values of e and z/ (see, 
for example, Fig. 0), relation ( |117|) warns that the corrections due to non-zero mass 
cannot be discarded. In particular, for the regime of very small energies and v = 
the term (p 2 /f 2 e/ln(l/e) and obviously exceeds the propagator correction 

F( 2 ) r-u g 2 . The conclusion is the following: in order to be within the accuracy set 
above one has to use the expression w k = y/k 2 + m 2 for the energy, though it is 
enough to use the residue J(p) of the instanton solution of the massless theory. 

For the residue of the instanton propagator and the propagator correction a similar 
analysis can be carried out. A general expression for the Fourier transform of the 
instanton propagator in the massive theory is of the form 

2 h(k 2 p 2 )q 2 p 2 ,sp 2 ,m 2 p 2 ) 

where s — (k + q) 2 . Expanding the function in the numerator in powers of m 2 p 2 , 
h(k 2 p 2 , q 2 p 2 , sp 2 , m 2 p 2 ) = h {k 2 p 2 , q 2 p 2 , sp 2 ) + (pm) 2 h 2 {k 2 p 2 , q 2 p 2 , sp 2 ) + ... 
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(logarithmic terms, like p 2 m 2 ln(p 2 m 2 ), may also appear), we obtain that the double 
on-mass-shell residue of the propagator can be written as 

P (m W 2 ) = p 2 h (OAs#p 2 )+p 2 (pm) 2 [h 2 (0,0,s # p 2 ) 
- h' 01 (0,0,s # p 2 )-h' 02 (0,0,s # p 2 ) 

where # = aa, ab, bb, s# is given by Eqs. (Ejj), (^Sj) and the upper index (m) refers to 
the non-zero mass case. It is easy to see that the residue of the instanton propagator 
in the massless theory is then given by 

P(sp 2 )=p 2 h (0,0,sp 2 ). 

A typical term in the propagator correction is 
dk dq 

where K(wi l Ti,w q x) is defined by Eq. (|116|) . By a simple analysis one gets 



/t/k da 
I cA K(w k T U w k x)P {m) (p 2 s#(K q))ir(w q r 2 , w^h 
W\c Wet 



/c/k c/q 
I cA K(w k T ll w i a)P (m \p 2 s#{k ) q))fT(u; q T 2 ,w q x)/ c 



P 2 



/c/k dc^ 
IK(w k Ti, w k x)P(p 2 s # (k, q))ir(u7 q r 2 , w q %) I 



+0(^(pm) 2 ,(pm) 



1 



r 4 



Again, from our numerical results we see that the corrections are small comparing 
to the propagator correction F^ 2 \e, u), calculated with / and P(p 2 s), and can be 
neglected without loss of accuracy. Similar to estimate (|117|) , if instead of w k = 
y/k 2 + m 2 and s#(k, q) the expressions |k| and s# (k, q), Eq. ([82l) , respectively, are 
used, then one can easily show that 



P 2 



IK{vw u w k x)P{p 2 s#{k ) q))^(w q T 2 , w^x) 1 

= P 2 J Ij^j/Kdkln, |k| X )P(p 2 ^ 0) (k, q))K(|q|r 2 , |q| X )/ 

+ {^)0 ((mr) 2 , (mp) 2 In ^ (mp) 2 ) , (118) 

where we assumed that r ~ T\ ~ r 2 . Since for a certain range of e and v the saddle 
point value f ~ 1, Eq. (|118|) gives an indication that the mass corrections can be 
comparable to the propagator correction ~ (p/f ) 4 ln(p 2 /f 2 ) and should be taken 
into account. 

Now we can formulate a procedure of calculation of the leading and next-to-leading 
order corrections to the function F(e, v) as follows: in expressions for W and 
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1) the residue of the instanton solution I(p), given by Eq. ( [79|) , and the residue of 

the instanton propagator P(p 2 s), given by Eq. fl87D , of the massless theory are 
substituted; 

2) the energy w k and the s- variables s aa , s a b and s^, given by Eqs. fl24|) , fl25l) (or by 

by Eqs. (jSjJ), Q8"5D) of the massive theory are used. 

We have shown that this procedure is consistent and gives results which are within 
the accuracy, set by Eq. floT|), provided p 2 (e, v) and p 2 (e, z/)/f 2 (e, v) are small in the 
range of our calculations. From the results of our numerical computations in Sect. [7| 
we checked that this condition holds. 



9 Discussion and conclusions 

In the present paper we have analyzed the multiparticle cross section of the shadow 
processes induced by instanton transitions in the simple scalar model with the action 
given by Eq. (|59]). We calculated the exact analytical expression for the on-shell 
residue of the propagator of quantum fluctuations in the instanton background. Using 
this result we calculated the propagator correction (i.e. the next-to-leading order 
correction) to the suppression factor Fj A of the multiparticle cross section. 

The leading and next-to-leading (propagator) corrections to the function F(e, u) 
were calculated semiclassically, by evaluation of the leading and next-to-leading terms 
in Eq. (0), respectively, at the saddle point. 

Within the accuracy of the approximation, considered in this article, the saddle 



point equations are derived from the leading terms, Eqs. (|T3|), Q14]), in the "effective 
action". It turned out, that for such equations the saddle point solutions for the 
dimensionless parameters (mp) 2 , (mr) 2 , {mx) 2 an d 7 in the integral (O) exist only 



for a certain region of e = E/E sph and v = N/N sph , shown in Fig. [| 

For general e and v from this region we computed the leading and propagator 
corrections to the function F(e, v) numerically. For v = we derived the analytical 
expressions for F^(e, 0) and F^(e, 0), Eqs. fllOID , Q - in terms of the saddle 
point solution f (e) of Eq. (|103|) . For the regime of very small energies, namely when 
relation (|105|) holds, and v = we obtained explicit expressions for the saddle point 



solutions and the leading and propagator corrections to F(e, 0). Recall that according 
to conjecture (|II]) this is the function F(e, 0) which characterizes the cross section of 
the instanton induced processes with a few initial particles. 

The corrections FW(e, u) and F^ 2 \e^ v) approach zero value when e — > 0. However, 
apparently there is no clear expansion parameter. At v — one can see that at 
energies for which f(e) 1 



F«(e,0) = 2v ^ +..., (119) 
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F^(e,0) = Fgl f] (e,0) + ... = -Aa 2 v^yn^ + .... 

Hence, the combination (p 2 /f 2 ) ln(p 2 /f 2 ) can play the role of the small expansion 
parameter. For very small energies from Eqs. ( |106| ), ( |107| ) we obtain that 



p 2 M) 

p\e,0) p~ 2 (e,0) 
r 4 (e) 11 f 2 (e) 



1 + 



lnln(l/e) 




We see that in this regime the expansion parameter is eJ ln(l / e) . With these relations, 
of course, results ( |108| ) and ( |111| ) are easily recovered from ( |119| ) (recall that v = 192n 2 



and c*2 = 3/2). 

The range of validity of the next-to-leading order of approximation of the function 
F(e, v) was estimated by comparing our results with numerical computations of the 
function F(e, v) in Ref. B for values of e and v for which the latter can be translated 
to the case of shadow processes. Such translation and the subsequent comparison 
of the results was done for the periodic instantons. The comparison shows that our 
perturbative results do not differ significantly from the exact ones for e < 0.25 or, 
equivalently, for v < z/ p (0.25) ~ 0.2. Thus, the intesection of the region e < 0.25, v < 
0.2 with the region for which the saddle point solutions exist (Fig. |3|) can be regarded 
as a rough estimate of the range of validity of the next-to-leading approximation. 

We would like to mention that, actually, in Ref. the function ^(e, v) (see 
Eq. ( |112| )) was calculated for points (e, v) such that 0.5 < e < 3, 0.2 < v < 1 and 
< Tie, v) < 0.6. For the range < e < 0.25 and < v < 0.2 and away from the 
line of periodic instantons methods of Ref. do not allow to obtain the value of T . 
Therefore, at the moment our perturbative calculations are the only ones which give 
quantitative behaviour of the supression factor in this range. 

We should stress that our result for the propagator correction is only the leading 
term of the expansion in powers of imp) 1 and p 2 /r 2 . This level of accuracy is de- 
termined by the fact that in formula flBTD for the action of the constrained instanton 
terms C(m 4 p 4 ) were not taken into account. In Sect. || we formulated the procedure 
of calculation of the constraint independent terms of the leading and next-to-leading 
corrections. The procedure essentially relies on the condition that the saddle point 
values for p 2 and p 2 /f 2 are small in the range of (e, v) under consideration. We checked 
that for the numerical saddle point solutions found the condition holds true. This 
conclusion is also verified analytically for the regime of very small energies, namely 
when relation (|105|) is valid. This shows the consistency of the procedure of calcu- 



lation. From this discussion it is clear that before calculating next-next-to-leading 
corrections to the function F(e, u) one should compute terms 0(p 4 ) in the leading 
and propagator corrections. 
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We also proved the cancellation of terms singular in the limit v — > in the 
propagator correction F^ 2 '(e, v) in a rather general context. The proof is essentially 
based on the general structure of the W^ l > and terms, the saddle point equations 
and the factorization property of the asymptotics of the propagator residue, following 
from a general formula of Ref. ||18|| . As we have explained, the problem of singularities 
in v is closely related to the problem of quasiclassical evaluation of contributions of 
initial states and initial- final states. In Sect. 14] we discussed this issue within the 



approach proposed by Mueller in Ref. [0. Namely, we calculated three leading 
terms of the asymptotics of the instanton propagator at large s and showed that with 
the appropriate choice of the propagator constraint the sins- and s-terms cancell out. 



According to Ref. |L2] , with such propagator the problem of semiclassical calculation 
of contributions due to initial states and initial-final states can be tackled properly. 

We expect that our results may provide some insight to the understanding of the 
structure of the suppression factor of the cross section of instanton induced processes 
in more realistic models, like QCD and the electroweak theory. They may also help to 
describe some features of the behaviour of such cross sections for of energies E <C E sp h, 
which are of interest for some planned high energy experiments [[UI . 
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Appendix A 



Evaluating Eq. (|49|) at the saddle point solution (^) - (|101|) with the residue of the 



instanton propagator given by expression (|88| ) we obtain the propagator correction to 
the function F(e, v) at v — > by straightforward calculation. Before presenting the 
result let us introduce some definitions and explain the origin of some terms. 
First we write the s- variables, defined by Eqs. (f24|), (p5|), as 



s# (k, q) = 2kql # [a # (k, q) - cos 9] , 

where # = aa,ab,bb, k and q stand for |k| and |q|, respectively, l aa = 1^ = — 1, 
l ao = 1, and 9 is the angle between k and q. While calculating the function S#(k, q), 
given by Eq. ([52]), one encounters the integrals 



M# tl (k,q) = f sin 9d9 ln(a # (fc, q) - cos 9) = a# In ° # + - + ln(ai - 1) - 2, 
jo ct^t — 1 
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M#,2(k,q) 



sin 9d9(a#(k, q) — cos 6) ln(a#(A;, q) — cos i 



4 + 1 ,„ °# + 1 



; ln^-^ + a # ln(a^-l)-a # . (Al) 

The logarithmic terms in the integrands above are due to the logarithmic terms in 
the function P(p 2 s) (see Eq. (j88[)). To get contributions to the propagator correction 
the functions M# 7 i(k,q) and M# t2 (k,q) are to be integrated over k and q (see Eq. 
(|5l|)). Hence, we define 

^#,i(ri,r 2 ) = [™ —^e-^e-^ikqf^M^q), 

JO Wk W q 

where i — 1,2. In particular, we will need the function N ab ^ 2 m the case when one of 
its arguments goes to zero. We define 

^2 



V(T) = -±N abi2 {T,T 1 )\ n _ Q . 



It can be shown that 



k 3 dk 



e- w * T M ab , 2 {k,oo) 



'0 w k 

where M ab ^ 2 (k, oo) is calculated from (|A1[) : 
M abt2 {k, oo) = — — In — — — In — — — \ 1 + 2 In — ] . (A2) 



2A; 2 



m 



Ik 2 



m 



k 



m 



The logarithmic terms in the propagator residue P(p 2 s) also give rise to the intrgrals 
of the form 

r°° dk ( k \" 



"^ T ln-. 

m 



They can be written as 



d r°° dk ( k 



-W k T 



d(5 Jo Wk \m 



d 

d[3 



1 ( 2 



- 



7T \VTLT J 



'(3 + 1 



K 



13/2 



VtVT , 



/3=n 



Let us introduce the function 

1 d 



mr dp 



Kp{mT)\p =1 + §{mT) 



mr 



C E + \n — 



where $(mr) is given by Eq 
Ai(mr) 
A 2 (mr) 
A 3 (mT) 



Finally, we define the functions 



— *'(mr) + — $(mr) + -n(r), 
2 mr 2 

+ -JV«,,i(r,r), 

2 1 
$'(mr)$'(mr) $(mr)$'(mr) - $(mr)f(mr) + -N bbt2 (T,r). 
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The partial contributions to the function F(e, v) in the limit v — > are equal to 

(A3) 



-32aivp 6 



(X 



71 P~ 2 
30 + n (x-r) 2 



_ 8«i e 3 
" v 2 $'(f) 



1 109 
2hi-r - 12$(f) +21n(-u$'(f)) -6C £ + 41n2- — 

32«!t;p 6 7 [ / p 2 71 , \ . 



(x - ^ ) 3 



(A4) 
(A5) 



16«i e 3 r 1 49 

In — — 8$(f ) + In (— u$(f)) — 4Ce — -r^ + 2 In 2 



lnr 

,(2) 



v 2 &(f) 
A{f) 



$'(f 

F ( i ;i /} (e, 0) = -8cw 
p 2 28 



(^)f(ln^ + ^-ln2 



p 2 28 



(A6) 
(A7) 



ln^ + ^-ln2l$ 2 (f)+^ 3 (f) 



— Aa 2 v p 



t 7 1 
f 2 + 15 



ln^r + -^-ln2 }<S> 2 (?)+A 2 (f) 



?3 



8«i e 
v 2 $'(f) 

$ 2 ^) 
(<&'(f)) 2 



-4*(f)-2lBf-2^-l-ln2 + ^ 



f 2 



-4$(f) - In y - 2C E - — - In 2 



u ($'(f)) 2 



ln2 + ^-2C E -2-4$(f)ln^ + ^ 
15 4 $ 2 (r) 



(A8) 
(A9) 



where, as before, we use the notation v = 1927T 2 . Note that these are the terms 
lnp/(x — f) in Eq. ( |A3| ) and Eq. ( TA5| ) that give rise to singular terms ln(l/V). The 
latter are shown explicitly in Eqs. (|A4|), 

For r C 1 the expressions above simplify considerably. The functions Ai(mr) 
have the properties 



Ai{mr) 

$'(t) 
Azimr) 



In 2 + C(mr), ^j" 1 ') = In 2 - 1 + 0(mr) 



$ 2 (mr) 



In 2 + i + 0{mr) 



($'H) 2 

Using these relations it is easy to calculate the expressions for the partial contributions 
(TAID - ( |Xf ) for r < 1. One gets 



(er)< 



f 2 



2(1 + — + C(f 3 ) J In ^ - ^ - 12 In 



1 12 



4 
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7 m 
(/-/) 



+ 2\nv-l2C E - 



+ 



109 



30 



v 



+ In v - 8C E - 



f 2 



+ 0(f") ln- 



^ r 2 



8 In 



r 
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4«i ,„„, 



4a 2 



30 

4 f 41' 

— - -41n- -AC E + — 
f 2 2 £ 30 



l4(81n^ + 8^-|) + ^) 



(A10) 

(All) 
(A12) 
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